Viscous evolution of the rapidity distribution of matter created in relativistic 

heavy-ion collisions* 
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Longitudinal hydrodynamic expansion of the fluid created in relativistic heavy-collisions is con- 
sidered taking into account shear viscosity. Both a on-vanishing viscosity and a soft equation of 
state make particle distributions in rapidity narrower. The presence of viscosity has dramatic con- 
sequences on the value of the initial energy density. The reduction of the longitudinal work and 
dissipative processes due to the shear viscosity, increase the total entropy and the particle multiplic- 
ity at central rapidities. The total energy in the collision, dominated by the longitudinal motion, is 
conserved. Viscous corrections make the longitudinal velocity of the fluid to stay close to the Bjorken 
scaling v z = z/t through the evolution. At the freeze-out viscous corrections are the strongest for 
non-central rapidities. 
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I. INTRODUCTION 

Properties of hot and dense strongly interacting matter 
can be studied in ultrarelativistic nuclear collisions. The 
modeling of the evolution of the dense collective phase 
is most commonly undertaken within the framework of 
relativistic fluid dynamics [3, H, Hj] . A clear evidence of a 
collective behavior of the system created in a collision is 
given by the observation of transverse elliptic flow of the 
produced particles. Collective flow arises naturally dur- 
ing a hydrodynamic evolution. In ultrarelativistic heavy- 
ion collisions the movement of the matter at the initial 
stage is dominated by the expansion in the longitudinal 
direction. Most of the experimental data at the Rela- 
tivistic Heavy- Ion Collider (RHIC) are restricted to the 
central rapidity region. Therefore, hydrodynamic mod- 
els often assume a simplified geometry of the fireball with 
a Bjorken boost invariant flow in the longitudinal direc- 
tion and concentrate on the dynamics in the transverse 
directions with azimuthal symmetry for central collisions 
[H HI or azimuthally asymmetric geometry for colli- 
sions of nuclei at non-zero impact parameters @, Hi & El • 
Only few calculations consider a fully three-dimensional 
evolution of the fluid [111, LL2J, |l3( . Results of the hydro- 
dynamic evolution are sensitive to the chosen equation of 
state. Values of the Hanbury Brown-Twist (HBT) radii 
observed at RHIC [1J| seem to exclude equations of state 
with a strong first order phase transition, or even with a 
more prominent soft point. The hydrodynamic evolution 
depends on the initial temperature and its profile, on the 
chosen equation of state, and on the freeze-out temper- 
ature. Particles emitted at the hydrodynamic freeze-out 
can still rescatter, which modifies somewhat their spectra 
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and elliptic flow [H M, M, [T^ . 

Instead of a fully three-dimensional calculation only 
longitudinal expansion of the matter created in a collision 
can be considered 0, [l8| . Recently the rapidity distribu- 
tions of pions and kaons produced in ^/s = 200AGeV 
Au-Au collisions measured by the BRAHMS Collabora- 
tion have been analyzed in a 1+1 dimensional hy- 
drodynamic model [20]. The final rapidity distribution 
of mesons has been found to be sensitive to the initial 
energy density distribution in rapidity and to the chosen 
equation of state. Experimental data indicate that the 
boost invariant Bjorken scaling solution is not realized at 
RHIC energies. Results of the 1+1 dimensional longitu- 
dinal fluid dynamics for the meson rapidity distributions 
favor a soft equation of state and require a Gaussian ini- 
tial energy density distribution in the space-time rapidity 

Analysis of the elliptic flow in the momentum distribu- 
tion of particles as a function of the initial eccentricity of 
the source points to the possibility of a noticeable effect 
of the shear viscosity of the fluid [U, [22j • Quantitative 
studies have to take into account many effects, initial 
conditions, the equation of state, the freeze-out temper- 
ature, possible final rescattering, therefore estimates of 
the viscosity coefficient are still under debate 0, [2~j |. 
Theoretical estimates of the ratio of the shear viscosity 
coefficient r\ to the entropy densitys range from a con- 
jectured lower bound r]/s = l/4n |25[ to r\ ~ s [2ol |27|. 

The role of the shear viscosity in the dynamics is most 
important during the early evolution of the system, when 
the velocity gradients are the largest. Gradients of the 
Bjorken flow give rise to corrections of the pressure ten- 
sor in the fluid. The transverse pressure increases, and 
the fluid expands faster in the transverse directions, this 
leads to stronger transverse flow and also to the satura- 
tion of the elliptic flow [H, H H, M, SH H, [H . All re- 
cent calculations using viscous relativistic hydrodynam- 
ics assume boost invariant Bjorken flow in the longitudi- 
nal direction and study the transverse development of the 
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fluid in azimuthally symmetric or asymmetric conditions. 
Longitudinal pressure is reduced, and hence so is the lon- 
gitudinal flow of the fluid. The fluid cools slower, at least 
until substantial transverse flow builds up. Weaker longi- 
tudinal expansion and entropy production due to dissipa- 
tive evolution require an adjustment of the initial entropy 
(temperature). As a result, the lifetimes of the plasma in 
the viscous and ideal fluid evolutions are similar. Finally, 
let us note that viscous corrections to the distribution 
functions at the freeze-out modify the spectra, and the 

HBT radii mil. 

Enhanced transverse pressure and the resulting modi- 
fication of the spectra at central rapidity, raise the ques- 
tion of possible modifications of the fluid dynamics in the 
longitudinal direction due to viscosity. Since the boost 
invariant scaling solution is not applicable at RHIC ener- 
gies, quantitative description of the energy flow and en- 
tropy production in the fireball should take into account 
a fully three dimensional geometry. Such a task within 
the dissipative fluid dynamics has not been accomplished 
yet. In the following we consider a simpler problem of 
the evolution of a longitudinally expanding non-boost in- 
variant fluid with viscosity. One expects that reduced 
longitudinal work in a viscous fluid generates narrower 
particle distributions in rapidity. The author is aware of 
only one, two-decades old work, considering this problem 
[HI, where first-order viscous hydrodynamics has been 
applied. At lower energies no strong effect of viscosity on 
the longitudinal expansion has been observed [35|. 

In the following we use BRAHMS data [19| to con- 
straint meson rapidity distributions after freeze-out. 
Solving viscous hydrodynamics in the 1+1 longitudinal 
geometry, we find a significant reduction of the initial en- 
ergy density when viscosity is taken into account. Also 
the longitudinal flow is modified, viscosity reduces the 
flow and counteracts the acceleration due rapidity gradi- 
ents of the pressure. At the Large Hadron Collider (LHC) 
energies we find that the modification of the longitudinal 
dynamics due to shear viscosity leads to an increase of 
the rapidity range where the Bjorken scaling flow applies. 



and the pressure are related by the equation of state. The 
energy density e(t, z) and the longitudinal velocity com- 
ponent v(t, z) are functions of time t and the beam axis 
coordinate z. Instead of the time and the z coordinate 
it is preferable to use the proper time r = \/t 2 — z 2 and 
the space-time rapidity 
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Hydrodynamic equations d^T^ — can be written as 

MM 



(e + p)Du" = Vp - A£ v Q 7r" a + tt^Duu (2.4) 



and 



De = -(e + p) V^u" + -n^ < V M u„ > 



(2.5) 



where 



D = = cosh(r - 9)d T + Smh(F 6 ^ d e , (2.6) 



< Vpii,, >= V^Uv + V„m m - -A^V a u a , (2.7) 



(- sinh YfC, 



-d x ,—d y , — coshYJC) (2. 



with A, 



u^u v and 



K = sinhfY - 6)d T + COSh(F 9) d e . (2.9) 

T 

The equations of the second order viscous hydrodynamics 
(|2.4p and (|2.5p are supplemented with a dynamic equa- 
tion for the stress tensor [H, H 13, H |M HJ 



T) < W > -27V7T 



(2-10) 



II. LONGITUDINAL HYDRODYNAMIC 
EQUATIONS WITH SHEAR VISCOSITY 

We consider a baryon-free fluid with non-zero shear 
viscosity. The energy-momentum tensor is the sum of 
the ideal fluid component and the shear tensor n 1 *" 



T» v = (e + p)u»u u - + 



(2.1) 



e and p are the local energy density and pressure of the 
fluid, and 



u» = 7(1,0,0,1;) 

= (cosh Y, 0,0, sinh Y) 



(2.2) 



is the four- velocity of the fluid element (7 = 1/ \/l — v 2 ) , 
and Y = i In fi^H is its rapidity. The energy density 



77 is the shear viscosity coefficient and tv is the relaxation 
time of the stress tensor, w'"' = A /ia A I//3 (d a up — dpu a ) 
is the vorticity of the fluid; it is zero for the longitudinal 
flow considered here. The relaxation time and the vis- 
cosity coefficient can be estimated from microscopic mod- 
els, considering equilibration processes [lH, [H, 13, 5l[ . In 
this work we take several values for the ratio of the viscos- 
ity coefficient to the entropy rj/s, and we drop viscosity 
effects for temperatures below 130MeV, the relaxation 
time is [26J t^/t] = 6/Ts, unless specified otherwise (T 
is the local temperature). For short relaxation times the 
stress tensor relaxes fast and stays close to the Navier- 
Stokes value 7r Miy = 77 < V^u" >. The dependence on 
the initial value of the stress tensor and on the relax- 
ation time is discussed in Section IVTl 

For a fluid expanding only in the longitudinal direction, 
with the energy density and the velocity constant in the 
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transverse plane, the stress tensor can be written using 
one scalar function II 



/ - sinh 2 Y 0- sinh Y cosh Y \ 
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n 



\-sinhycoshy -cosh 2 F / 



(2.11) 



and viscous hydrodynamic equations take the form 



(e+p)DY = -tCp + UDY + KH 
De = (e + p)KY - UKY 



DU 



{-■qKY - nj/r, 



(2.12) 



The above equations are solved numerically in the t-6 
plane, starting from some energy density e(ro,9) at the 
initial proper time of the evolution To = l.Ofm/c. For the 
initial fluid rapidity we always take the Bjorken flow 



Y(t ,6) = 9. 
Assuming boost invariance, i.e. 

Y(t,6)=6, e(r,0)=e(r) , p(r,0) = 
hydrodynamic equations simplify to [37t [3 

e + p- n 



(2.13) 
p(r) (2.14) 



de 
dr 
dU 



(4t7)/(3t) - n 



(2.15) 



III. EQUATION OF STATE 

The equation of state determines the evolution of the 
fireball. In the following we take a parameterization 
of the equation of state proposed by Chojnacki and 
Florkowski [42] . It is an interpolation of the lattice data 
at temperatures above T c — 170MeV and of an equation 
of state of noninteracting hadrons at lower temperatures. 
The two limiting formulas are joined smoothly with only 
a slight softening of the equation of state around the 
critical temperature (Fig. [I]). This minimally softened, 
realistic equation of state has shown itself suitable for the 
description of the transverse expansion of the fluid and 
of the buil d up of the elliptic flow and gives reasonable 
HBT radii [1||43|. 

In this section we consider the hydrodynamic longitu- 
dinal expansion of an ideal fluid. Equations for the fluid 
rapidity and energy density are obtained from Eqs. (|2.12p 
setting II = 0. This problem has been discussed recently 
in the context of heavy-ion collisions at RHIC [2(| ■ For 
completeness we study the effect of the equation of state 
we use on the longitudinal expansion. For that purpose 
we compare the evolution using the equation of state of 
Chojnacki and Florkowski with an evolution based on a 
relativistic gas equation of state p = ie. 
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FIG. 1: Square of the velocity of sound as a function of the 
temperature for an equation of state interpolating between 
the hadron gas and the quark-gluon plasma expressions [421 ]. 



At the initial time To = l.Ofm/c the energy density is 

e(T ,6) = e exp(-9 2 /(2a 2 )) . (3.1) 

The initial energy density eo and the width of the ini- 
tial rapidity distribution a are parameters adjusted to 
reproduce final meson rapidity distributions. The freeze- 
out takes place at the temperature Tf = 165MeV, this 
high freeze-out temperature is the same as the chemical 
freeze-out temperature [3, EH . The spectrum of parti- 
cles emitted with four-momentum = [E, q) is given by 
the Cooper-Frye formula [2(| EH 



E 



d 3 N d 3 N 



1 



d 3 q dyd 2 q± (2tt) 3 



dS^ftfuJ , (3.2) 



f(E) — e~ E l T s is the thermal distribution (Boltzmann 
distribution for simplicity), and y denotes the rapid- 
ity of the emitted particle. The element of the hy- 
persurface of constant freeze-out temperature is ciE^ = 
S (t\0) sinh 6+t(9) cosh 9, 0,0, r (9) cosh 9-t{6) sinh 9), 
S = kR\ u is the transverse area of the fireball in central 
collisions and t(9) is the line of constant temperature Tf 
in the t-9 plane. Particle distributions in rapidity are 
obtained by integration over the transverse momenta q± 
inEq. {321 



dN S 
dy 4-7T 2 

(2r< + 2£ 2 + to 2 )£ 

mcosh(y — Y;{9)) 



(t(9) cosh(y - 9) - t (9) sinh(y - 9)) 



exp 



T 



f 



d9 



(3.3) 



m is the meson mass, YA9) = Y(r(9),9) is the fluid 
rapidity at the freeze-out hypersurface, and 



cosh( 2/ - Y f {6)) 



(3.4) 
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FIG. 2: (color online) Rapidity distribution of mesons calcu- 
lated using a realistic equation of state (solid line) and using 
a relativistic gas equation of state for three different initial 
conditions (dashed, dotted and dashed-dotted lines), Table I. 
Data are from the BRAHMS Collaboration fijl. 



The above expression neglects the transverse expan- 
sion of the fluid at the freeze-out. This has only a small 
effect on rapidity distributions, the distribution should 
be narrower. Pions and kaons come to a large extend 
from secondary decays of resonances. The emission takes 
place in two stages, first an emission of a heavy resonance 
according to Eq. (|3.3[) and then the decay of the reso- 
nance into pions (kaons). The emission of resonances and 
their decay is also influenced by their transverse expan- 
sion. A hint on the spread in rapidity of the decay prod- 
ucts of resonances is given by charge balance correlations 
pTl . I4H [49j . Narrow charge balance functions indicate 
that decay products of a resonance are only 0.5 unit of ra- 
pidity away from the parent resonance, convoluting this 
distribution with the spread in rapidity of the emitted 
resonances one obtains a distribution of half- width simi- 
lar as for the emission of direct pions in Eq. (|3.3p . Since 
75% of pions come from resonances at Tf = 165MeV [50] |. 
we multiply the distribution from (|3 . 3|) by a factor 4 to 
account for all pions, direct and from resonance decays 
(the factor is 1.7 for kaons). 

The parameters eo an d o have been adjusted for the 
calculation using a realistic sound velocity (Fig. [1} to 
reproduce the width and normalization of the observed 
pion distribution. The resulting meson distributions are 
similar as the ones observed experimentally. When using 
a relativistic gas equation of state (p = ^e) one always 
gets a meson distribution in rapidity that is too wide. We 
present three calculations with different initial widths a 
and with the initial energy densities £q adjusted to re- 





eo 

GeV/fm 3 


a 


r(0) 
fm/c 


EOS [42] 


16.9 


1.05 


14.8 


p = |e case I 


71.5 


1.05 


14.8 


p = |e case II 


102 


0.8 


15.7 


p = |e case III 


50.8 


1.5 


13.7 



TABLE I: Parameters of the initial energy density distribu- 
tion (|3.1[) . for ideal fluid calculations, one using a realistic 
equation of state and three calculations using a relativistic 
gas equation of state. In the last column is shown the lifetime 
of the system until freeze-out. 
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FIG. 3: (color online) Difference between the flow rapidity 
of the fluid and the Bjorken value, calculated using a realis- 
tic equation of state (solid line) and using a relativistic gas 
equation of state for three different initial conditions (dashed, 
dotted and dashed-dotted lines), Table [I] 



produce dN / dy for central rapidities only (Fig. [2]) . The 
final meson distribution is much wider than the initial 
energy density distribution in all cases. It is due to the 
breaking of the Bjorken scaling of the longitudinal flow, 
Y(t 7 9) > 6 (Fig. [3]). Fast moving fluid elements emit 
mesons in the far forward and backward rapidities. Con- 
cluding this section, we confirm the findings of Ref . [2(| ■ 
A hard equation of state never works; a narrow initial 
distribution of the energy density in rapidity leads to a 
strong acceleration of the longitudinal flow, wider initial 
distributions are incompatible with the narrow final me- 
son distributions. Only by imposing a softened equation 
of state, experimental pion and kaon distributions can be 
approximately reproduced. In the following we take the 
realistic equation of state from Ref. [42| and study the 
effect of non-negligible shear viscosity. 
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FIG. 5: Relative viscous corrections TL/(e+p) at the freeze-out 
as a function of space-time rapidity (Eq. 14. 4p . 



FIG. 4: (color online) Rapidity distribution of mesons calcu- 
lated using a realistic equation of state and viscosity 77/s = 0.2 
(solid line) and using ideal fluid hydrodynamics (dashed line). 
The dotted line denotes the results of a viscous hydrodynamic 
evolution, but neglecting the viscous corrections to the par- 
ticle emission at freeze-out (Eq. 14.3) 1 . Data are from the 
BRAHMS Collaboration [H. 



IV. DISSIPATIVE LONGITUDINAL 
EXPANSION 

Following the hydrodynamic evolution with viscosity 
requires the solution of the coupled equations (|2.12p , for 
the fluid rapidity Y, the viscous correction II, and the 
energy density e. The initial conditions are given in Eqs. 
p.l[) and (|2 . 13|) . For the initial viscous corrections we 
take U(to,8) = p(tq,8). This means that at the initial 
time the effective pressure is maximally anisotropic in the 
fluid rest frame, with zero longitudinal pressure. We vary 
the value of the shear viscosity coefficient 77/s, and in each 
case we adjust the parameters of the initial energy density 
distribution (|3.1[) to reproduce BRAHMS data for the 
pion distribution in rapidity. Viscous corrections modify 
the momentum distribution functions in the fluid. We 
assume the same modification for all the particle species 
in the fluid 



f(q)+Sf(q) = f(q)[l + 



2T 2 (e+p) 



(4.1) 



The correction to the particle distribution function is 
proportional to the ratio of the viscous correction n to 
the enthalpy e + p. Viscous corrections to the distribu- 
tion functions at the freeze-out modify the Cooper-Frye 
formula 



dy 



dN_ 

dy 



dSN 
dy 



(4.2) 



To the expression 



one has to add 



dSN 
dy 



S 
4tt 2 



(t{9) cosh(y -6)-t (9) sinh(y - 9)j 



[12£ 5 + 5£ 3 m 2 



12£ 4 m + £ 2 m 3 - sinh(y - Y f (9)) 
(24£ 5 -I- 12£ 3 m 2 + 24£ 4 m + 4£ 2 ™ 3 + £m 4 )] 
n ( mcosh(y-Y f {6)) 



2T 2 {e+p) 



exp 



T 



f 



d9 



(4.3) 



In Fig. [4]is shown the result of this procedure for the case 
77/s = 0.2. Again the direct meson spectra are multiplied 
by a factor 4 for pions and 1.7 for kaons, to account 
for the expected ratio of all mesons to directly produced 
mesons [50j |, at the chosen freeze-out temperature. Pion 
emission at the end of the viscous hydrodynamic evolu- 
tion (solid line) is similar as observed experimentally. In 
Fig. 2] is also shown the meson distribution obtained us- 
ing the equilibrium distribution at freeze-out (Eq. 13. 3p 
(dotted line). Deviations from the full result (Eq. 14. 2[) 
is only noticeable at rapidities 3 unit away from central 
rapidity. It can be understood as due to an earlier freeze- 
out at large rapidities, which makes the relative viscous 
corrections 



H(t(0),0) 



e{t{6),0)+p{t{6),0) 



(4.4) 



at freeze-out more important (Fig. [5J. Accordingly one 
expects the largest correction from viscosity to particle 
spectra, elliptic flow, and HBT radii at large rapidities. 

We have noticed (Section IIII[) that longitudinal gra- 
dients of the pressure cause the acceleration of the mo- 
tion of the fluid in the beam direction. On the other 
hand, viscosity reduces the longitudinal motion of the 
fluid, and as a consequence reduces its expansion in the 
longitudinal direction. Shear viscosity prevents large gra- 
dients of the velocity field to develop. One can compare 
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FIG. 6: (color online) Difference between the flow rapidity of 
the fluid and the Bjorken value, calculated for an evolution 
with shear viscosity coefficient rj/s = 0.2 (solid lines), for an 
ideal fluid with a realistic equation of state (dashed lines) 
and using a relativistic gas equation of state (case I) (dashed- 
dotted lines). 



FIG. 7: (color online) Difference between the flow rapidity of 
the fluid and the Bjorken value at the freeze-out hypersur- 
face, calculated for an evolution with shear viscosity coeffi- 
cient rj/s = 0.2 (solid line), for an ideal fluid with a realistic 
equation of state (dashed line) and using a relativistic gas 
equation of state (case I) (dashed-dotted line). 



■q/s 


eo 


a 


r(0) 




GeV/fm 3 




fm/c 





16.9 


1.05 


14.8 


0.1 


9.8 


1.18 


14.1 


0.2 


5.6 


1.8 


13.1 


0.3 


4.0 


1.86 


12.4 



TABLE II: Parameters of the initial energy density distribu- 
tion (|3.1|) for hydrodynamic calculations with several value of 
the shear viscosity coefficient. In the last column is show the 
lifetime of the system until freeze-out. 



the rapidity of the fluid with and without shear viscos- 
ity (Fig. [5]). For rj/s = 0.2, the flow stays close to the 
Bjorken one during a substantial part of the evolution. 
At the freeze-out hypersurface the flow is still Bjorken- 
like for \0\ < 1.8 for rj/s = 0.2 (Fig. [7]). The acceleration 
from pressure gradients and the deceleration from viscos- 
ity approximately cancel for this choice of parameters. 
Reduced longitudinal expansion with viscosity requires 
smaller initial energy densities to reproduce the final me- 
son distributions. In Table ITI1 are listed the initial energy 
densities and widths of rapidity distributions adjusted 
to reproduce the observed meson distributions for sev- 
eral values of the shear viscosity coefficient. The shape 
of the initial energy density is extremely sensitive to the 
dynamics of the longitudinal expansion (Fig. [5]). Obvi- 
ously the value of the initial energy density, or in other 
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FIG. 8: (color online) Initial energy density distribution for 
the ideal fluid hydrodynamic evolution with a realistic equa- 
tion of state (dashed line), for viscous hydrodynamic evolu- 
tions (solid lines), and for a relativistic gas equation of state 
(case I) (dashed-dotted line). 
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words the cooling rate from the longitudinal motion, is 
an important ingredient in the modeling of the transverse 
expansion of the fluid. Let us also note, that when re- 
duced initial energy densities are imposed, the lifetime of 
the system until freeze-out is only weakly dependent on 
the viscosity coefficient. 



V. COOLING AND ENTROPY PRODUCTION 

It is instructive to compare the cooling rate in our so- 
lution and in the boost invariant scaling solution. Finite 
extension in space-time rapidity makes the cooling rate 
faster. On the other hand, reduced velocity of sound and 
shear viscosity reduce the longitudinal work of the pres- 
sure, and slow down the cooling. In Fig. [9] is compared 
the cooling of the central region of a finite system to the 
cooling in the boost invariant case (Eq. I2.I5|) . To com- 
pare the solution in the 1+1 dimensional system to the 
boost-invariant one, the initial temperatures are fixed so 
as to give the same lifetimes of the two systems until 
freeze-out. Boost invariant solutions underestimate the 
values of the initial temperature (energy density) and of 
the cooling rate for an ideal fluid evolution. In a finite 
system additional cooling appears and the longitudinal 
flow is stronger than the Bjorken scaling flow. As the 
shear viscosity coefficient increases the velocity gradients 
in the dynamics are more and more constrained. The 
flow accelerates less. At rj/s = 0.2 the effects of the 
viscosity and space-time rapidity gradients of the pres- 
sure counterbalance each other and the flow is approxi- 
mately Bjorken-like. With stronger viscosity rj/s = 0.3, 
the Bjorken flow is decelerated and the cooling is slower 
than for the boost-invariant solution. It means that the 
gradients of the viscous correction n are larger than the 
gradients of the pressures itself and applicability of sec- 
ond order viscous hydrodynamics is questionable. 

Dissipative hydrodynamics conserves the total energy 
but produces entropy. The expression for the total energy 
of the system at proper time r is 

/oo 
d6 [e cosh Y cosh(F - 6) 
-OO 

+ O-n)sinhYsinh(Y-0)] • (5- 1 ) 

It is mainly composed of the kinetic energy of the longi- 
tudinal motion of the fluid. Small changes of the energy 
density at large rapidities can cause significant changes 
of the energy density for central rapidities. Softening of 
the equation of state and non-zero shear viscosity modify 
the dynamics at large rapidities, which leads to less cool- 
ing at 8 = 0, while the global energy of the fireball is un- 
changed. Dissipative processes drivin g th e system lo cally 
towards equilibrium produce entropy [39|, Eo, El, [52|, [53| . 
The total entropy of the system 
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FIG. 9: (color online) Time dependence of the temperature 
at the center for a longitudinally expanding ideal fluid fire- 
ball (dashed line) compared to the boost-invariant solution 
(dashed-dotted line). Same for the evolution with viscosities 
r\l s = 0.1, 0.2, 0.3 (solid lines, increasing r\ from top to bot- 
tom), and for the boost-invariant case with viscosity (dotted 
lines). For rj/s = 0.2 the solid and doted curves lie on top of 
each other. 




S(t) = tS / d9s cosh(Y 



(5.2) 



[fm/c] 



FIG. 10: (color online) Relative entropy production in the 
viscous hydrodynamic evolution (Eq. 15. 2[) (solid line), of the 
entropy density at central space-time rapidity (dashed line), 
and of the entropy from the boost-invariant Bjorken solution 
(Eq. 15. 3p (dashed-dotted line). All calculation with rj/s = 
0.2. 
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FIG. 11: (color online) Ratio of the dynamical shear viscosity 
correction II(t) to the Navier-Stokes value ILys(t) at 77/s = 
0.2 for two initial condition for the viscous corrections n(ro) = 
p(tq) (solid line) and for II(ro) = (Eq. 16. ip (dashed-dotted 
line), and for a reduced relaxation time r n — 1.5ri/Ts (dashed 
line) . 



FIG. 12: (color online) Relative entropy increase in the dis- 
sipative evolution at rj/s — 0.2 for two initial condition for 
the viscous corrections II(to) = p(to) (solid line) and for 
II(ro) = (Eq. 16. ip (dashed-dotted line), and for a reduced 
relaxation time t w = 1.5r//Ts (dashed line). 



increases with time, if shear viscosity is active (Fig. ITU1) . 
The increase of the entropy at central rapidity ts(t, 6) is 
responsible for an increase of the particle multiplicity in 
the central rapidity region. Estimates of the entropy pro- 
duction [lU based on the boost-invariant solution (Eq. 

EH) 



d(rs) 
dr 



n 

T 



(5.3) 



are very close to the 1+1 dimensional dynamical result, 
it happens if the shear viscosity is strong enough to con- 
serve the Bjorken flow in the evolution. 



VI. ROLE OF THE INITIAL STRESS TENSOR 

Second order viscous hydrodynamics introduces a dy- 
namical equation for the stress tensor. A crucial param- 
eter is given by the relaxation time r n . For small re- 
laxation times the viscosity correction stays close to the 
Navier-Stokes value 



IWt) = g 



(6.1) 



Such a behavior of viscous corrections has been con- 
firmed in numerical simulations of the transverse ex- 
pansion in viscous hydrodynamics with small relaxation 
times [23], HH , unlike in Ref . [5l| where a large value of 



the relaxation time was postulated. We check the effect 
of the initial value of II(t) on the evolution by compar- 
ing two scenarios, II(to) = p(tq) which corresponds to an 
anisotropic momentum (pressure) at the initial time and 
II(to) = which corresponds to initially locally equili- 
brated distributions. With the choice of the the relax- 
ation time = 6r)/Ts the two initial conditions lead to 
different results. As before, for each choice of the initial 
conditions and parameters we readjust the initial energy 
density distribution (|3.1|) to reproduce the final pion ra- 
pidity distribution. 

In Fig. [IT] is shown the ratio of the dynamical value 
of the viscous correction II(t) to the Navier-Stokes value 
(|6.1|) . After several fm/c the stress tensor, that is set 
initially to zero (dashed-dotted line), relaxes to and over- 
shoots the steady flow value (|6.1[) . For an initial stress 
tensor corresponding to the anisotropic effective pressure 
(solid line) the viscous correction overshoots the Navier- 
Stokes value almost immediately. Since the dissipative 
effects are the strongest at the early stages, the integrated 
entropy production is smaller for an evolution with ini- 
tial zero stress tensor than for an evolution starting with 
nonzero stress tensor (Fig. IT2|) . Even though, the ef- 
fect of the change of the initial viscosity correction can 
be counterbalanced by a suitable change in the initial 
energy density distribution, the amount of the entropy 
produced in the dynamics depends on the initial value 
of the dissipative corrections. We have also performed a 
calculation with a smaller relaxation time (dashed lines 
in Fig. fTTI and [T^)l . In that case the initial value of II(t) 
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FIG. 13: (color online) Rapidity distribution of pions (upper 
curves) and kaons (lower curves, x3) for LHC plateau- like 
initial conditions (Eq. I7.1|l . calculated for an evolution with 
shear viscosity coefficient rj/s — 0.2 (solid lines) and for an 
ideal fluid with a realistic equation of state (dashed lines). 
The dotted line denotes the result of a viscous hydrodynamic 
evolution, but neglecting the viscous corrections to the parti- 
cle emission at freeze-out (Eq. 14. 3[) . 

is less important, as it rapidly relaxes to the the steady- 
flow value IIjvs(t). The integrated dissipative effects are 
smaller. 



VII. EXPECTATIONS FOR LHC 

In this section we present some simple estimates of the 
effects of the shear viscosity on the longitudinal expan- 
sion at LHC energies. To get a rough estimate we set 
arbitrary the multiplicity of pions for central rapidity at 
twice the value observed for central collisions at RHIC. 
The initial energy density distribution in space-time ra- 
pidity is modified, it includes a plateau of width a v (l3j 

e(T ,9) = e exp(~{e~a p ) 2 e(\e\-<T p )/2cT 2 ) . (7.1) 

The width of the plateau is a v — 3.3 and takes all the 
increase of the rapidity range when going from RHIC to 
LHC energies. The parameter a is the same as at RHIC 
energies, and the energy density eo is adjusted to repro- 
duce the assumed final pion density ^| y= o- The meson 
distributions are shown in Fig. Q2] For the viscous evolu- 
tion the plateau in the final meson distributions survives. 
For an ideal fluid one gets eo = 16.9GeV/fm 3 , while for 
a shear viscosity rj/s = 0.2 one needs eo = 12.4GeV/fm 3 . 
The difference between the initial energy densities in the 



FIG. 14: (color online) Difference between the flow rapidity of 
the fluid and the Bjorken value at t = lOfm/c, calculated for 
an evolution with shear viscosity coefficient rj/s — 0.2 (solid 
line) and for an ideal fluid with a realistic equation of state 
(dashed line) for LHC plateau-like initial conditions. 



viscous and ideal fluid evolutions is not as big as for RHIC 
energies. At LHC energies, both the ideal and viscous 
fluid evolutions have a Bjorken scaling form in several 
units of central rapidity (Fig [T4"|) . The lack of space- 
time rapidity gradients in the distributions makes the 
evolution and cooling last longer, around 20fm/c. Even 
at the freeze-out the flow is Bjorken-like at central ra- 
pidities (Fig. [T5]). This observation justifies the use at 
LHC energies of thermal and h ydr odynamic models as- 
suming boost invariance [H, l54l l55j . The hydrodynamic 
evolution can be restricted to the 1+2 dimensional boost 
invariant geometry to describe matter in the 4-5 central 
units of rapidity. At the freeze-out hypersurface dissi- 
pative corrections to the momentum distributions have 
almost disappeared (compare the solid and dotted lines 
in Fig. HU). 

Since the cooling of the fluid due to the longitudinal 
expansion is slow (like in the Bjorken solution), a realistic 
modeling of the time scales and of the freeze-out hyper- 
surface must take into account the transverse expansion. 
The speed of the transverse expansion would determine 
the life-time of the system, while at RHIC energies longi- 
tudinal expansion (in space-time rapidity) is also impor- 
tant. We have also performed calculations using Gaus- 
sian initial energy distributions (|3.1[) with rescaled width 
parameters a for the increased LHC rapidity range. For 
a viscous fluid we find a broad region of rapidities where 
the Bjorken scaling flow survives through the evolution, 
\8\ < 2.5. For the ideal fluid the scaling is broken, but to 
a significantly lesser extend than at RHIC energies. 
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FIG. 15: (color online) Difference between the flow rapidity 
of the fluid and the Bjorken value at the freeze-out hyper- 
surface, calculated for an evolution with shear viscosity co- 
efficient rj/s — 0.2 (solid line) and for an ideal fluid with a 
realistic equation of state (dashed line) for LHC plateau-like 
initial conditions. 

VIII. CONCLUSIONS 

The evolution of a fireball of dense and hot matter cre- 
ated in heavy-ion collisions can be modelled as a hydro- 
dynamic expansion of a viscous fluid. We analyze the ef- 
fects of the shear viscosity on the longitudinal expansion 
of the matter. We solve numerically coupled evolution 
equations for the longitudinal flow, the energy density 
and the viscous corrections in a 1+1 dimensional geom- 
etry, corresponding to a rapid expansion in the beam di- 
rection. As a function of the space-time rapidity the dis- 
tribution of matter evolves slowly with proper time. The 
average density drops and the distribution gets wider. 
The last phenomenon takes place when the flow of the 
fluid gets stronger than the Bjorken one. At the freeze- 
out temperature the hydrodynamic stage finishes and 
particles are emitted thermally according to the Cooper- 
Frye formula. Experimental measurements of the distri- 
bution of mesons in rapidity [l!^ ] constrain the allowed 
distribution of the longitudinal velocities of the fluid el- 
ements. Correlation between space-time and momentum 



rapidities of the fluid, means that the space-time rapid- 
ity extension of the fluid must be limited and that its 
longitudinal flow cannot deviate significantly from the 
Bjorken flow. Shear viscosity counteracts the gradients 
of the velocity field. As a consequence it slows down 
the longitudinal expansion. At the freeze-out the energy 
density distribution in space-time rapidity is narrower 
and the longitudinal flow gets less accelerated than for 
the ideal fluid hydrodynamics. Fitting the initial energy 
density distribution to reproduce the final meson distri- 
butions, one observes a striking effect. With increas- 
ing shear viscosity coefficient the initial energy density 
of the fireball decreases significantly, from 16.9GeV/fm 3 
for an ideal fluid to 5.6GeV/fm 3 for r//s = 0.2. In the 
1+1 dimensional longitudinal geometry, this energy den- 
sity corresponds to an average over the transverse plane. 
Nevertheless, estimates of the maximal energy density 
reached in heavy-ion collisions at RHIC energies must be 
strongly revised down if shear viscosity is effective during 
the expansion of the fireball. This dramatic reduction of 
the initial density should also be taken into account in 
hydrodynamic models dealing with transverse expansion 
only, both in 1+1 and 1+2 dimensions. 

Depending on the balance of the acceleration of the 
flow from pressure gradients and deceleration from vis- 
cosity, the flow gets faster or slower than the Bjorken 
one. For some values of the parameters, effects of the 
shear viscosity and pressure gradients on the longitudi- 
nal flow of the fluid cancel, i.e. the flow stays close to 
the Bjorken flow. This could be an argument justify- 
ing models, which combine transverse viscous expansion 
with a Bjorken flow in the beam direction. When the ini- 
tial conditions are adjusted to reproduce the final meson 
distributions, we find that the freeze-out hypersurfaces 
are very similar, irrespective of the value of the shear 
viscosity coefficient, obviously the lifetime of the system 
is not sensitive to viscous effects cither (Table [TTJ) . At 
the freeze-out the viscous corrections (from the longitu- 
dinal flow) to the thermal distributions are small, except 
at large space-time rapidities. At LHC energies a sub- 
stantial rapidity plateau, where Bjorken scaling applies, 
is expected to appear. Shear viscosity helps to preserve 
it in a wider rapidity interval through the evolution. 
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